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Abstract. - Within standard treatments of the interacting electron-phonon system the elec- 
tron density of states (EDOS) shows no sign of the phonons when the system is in the normal 
state. On the other hand, emergence of a fine energy scale in the superconducting state makes 
the phonons readily observable in the EDOS and other spectroscopic quantities. Here we wish 
to re-examine the conditions under which phonon structure might be present in the normal 
state. Accounting for a finite, i.e. non-infinite, electronic bandwidth results in easily distin- 
guishable phonon structure in the EDOS. We argue that for narrow band metals, the resolution 
of photoemission spectroscopy is sufficient for resolving these structures. Manifestations of fi- 
nite band effects in KiCao are discussed. 



Introduction. - The observation of phonon structure in tunneling data for superconduct- 
ing materials became a very important diagnostic tool in the 1960's [1]. Following this initial 
success attempts were made (with some degree of success) to extract similar information in 
the normal state [2-4] . Part of the difficulty at the time no doubt was due to the fact that 
calculations using many-body Green functions predicted that no phonon structure should be 
present in the Electronic Density Of States (EDOS) [5,6]. When the bare EDOS is infinitely 
wide and featureless and when particle-hole symmetry is present, the dressed EDOS retains 
its constant bare band value and the phonon signatures drop out. On the other hand, as was 
studied extensively in the 1980's (see for instance [7,8]), phonon renormalizations do enter 
when the EDOS has significant energy dependence near the Fermi surface on the scale of 
phonon energies. 

More recently [9-14] it was realized that phonon structure can also arise when a finite 
band cutoff is applied in the EDOS on the scale of the bandwidth. In this letter we want 
to accomplish two goals: first we want to understand under what conditions these structures 
can be resolved in Angle Integrated Photo Emission Spectroscopy (AIPES) and/or tunneling 
even for relatively broad bands, and secondly, we wish to establish the scale of the structures 
for parameters more specific to a real narrow band material, in this case KaCgo- While our 
primary interest here is in the EDOS we also provide comparison with other spectroscopic 
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quantities such as the real and imaginary parts of the quasiparticle self energy and its optical 
counterpart, the memory function, both of which have phonon signatures even when the band 
is assumed to be infinite. It is well known that even when these structures are not directly 
visible, as in the optical scattering rate, they can be emphasized either by differentiation [15], 
or extracted by other more sophisticated "inversion" procedures [15,16]. In contrast, the 
phonon structures due to a finite bandwidth are more readily detectable. 

To see sizeable features in the EDOS, the electron-phonon interaction must be sufficiently 
strong. To be specific we consider the Eliashberg a^F{ijj) function appropriate to the in- 
termediate coupling strength superconductor KsCgo. We start with a large bare bandwidth 
W = 2.5 eV and find that for the clean system and at low temperature the phonon structure 
in the renormalized EDOS is very prominent in the range of phonon energies (keep in mind 
that for an infinite bandwidth the structure is absent). Furthermore, the presence of impurity 
scattering does not obliterate the structure; it merely reduces its size [12]. Reducing the bare 
bandwidth to Vl^ = 0.5 eV, a value more relevant for KaCeo, we find that all dominant features 
of the renormalized EDOS are determined by the electron-phonon interaction. 

Bare EDOS, the self energy and the renormalized EDOS. - We use the following model 
for the bare electronic band: Nq{£^) — iVgO {W/2 — ]^]) , where W is the bare band width 
and Q{x) is the step function. The constant iVo is fixed by normalization: Nq ^ 1/W. In 
this letter we retain particle-hole symmetry for simplicity, with the chemical potential at the 
center of the band, ^ = 0. 

The electronic self energy T,{z) — + 122(2;) is calculated from the Migdal equations 

formulated in the mixed real-imaginary axis representation [12, 17]: 
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where cOm = 7rT(2m — 1), to e Z are the fermionic Matsubara frequencies, and /(w) and n{Lu) 
are the Fermi and Bose distribution functions respectively. The electron-phonon interaction is 
specified in terms of the electron-phonon spectral function a^F(uj) (the Eliashberg function). 
The parameter F, which has the meaning of an impurity scattering rate, specifies the strength 
of the interaction with impurities. Note that finding E(z) requires a self consistent solution 
of eqs. Q-lO; which is a consequence of accounting for the finite bare band width in eq. 

The variable z can assume arbitrary complex values. At first, the solutions for 'S{z) are 
sought on the imaginary axis, at z = iLOm, where eq. ^ is simpler. Then, the function r]{iujm) 
is used to set up an iterative procedure to find S(ijj) just above the real axis (we use uj as 
shorthand for a; + iO+). While we will focus on the EDOS, which can be determined in AIPES 
experiments, the real axis self energy, both its real and imaginary parts, can be obtained from 
the angular-resolved photoemission. The accuracy of this technique has increased dramatically 
in recent years and properties of both new and traditional materials have been scrutinized 
(see, for example, refs. [18-21] for a flavour of the recent developments). 

Motivated by the electron-phonon interaction in the fulleride compound K3C60, we use a 
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three frequency model for the electron-phonon spectral function: 



3 , 3 



a' 
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with h = 0.3, h = 0-2, h = 0.5 and ^1:^2:^3 = 0.04 : 0.09 : 0.19 eV [22]. The interaction 
strength a is defined as the area under the a'^F{u!) curve. The mass enhancement parameter 
A is given by eq. l|2Jl with z = 0. For the forthcoming discussion we set A = 0.71. Then this 
model has a = 43.8 meV and c^in = 102.5 mcV, where cui^ is the logarithmic frequency [6], 
a convenient parameter to quantify the phonon energy scale. In the last section we consider 
different values of A as well. Finally, below we use W — 2.5 eV. 

The renormalized density of electronic states (or density of states for quasiparticles) is 
defined by 



where A(^,lu) = —ImGretiS,,'-^)/''^ is the electronic spectral density and Greti^,^) is the 
retarded Green function. It can be expressed in terms of the function rj = rji + irj2 of Eq. 
(O through N{uj)/Nq{Q) — —ri2{uj)/iT. In AIPES the measured intensity, /(w), is the product 
of the renormalized EDOS and a thermal factor f{Lo) which provides a cutoff at T = 0, 
convoluted with the instrument resolution R{ijj) which can be taken as a Gaussian with half- 
width of 2.8 meV [20]. With particle-hole symmetry it is convenient to define a symmetrized 
quantity /s(w) = [/(w) + I{—uj)\ /2 from which f{uj) drops out [23], leading to 



In fig. n we show results (solid curves) for Ig vs lo and compare with the density of states 
N{u}) (dotted curves). The main frame is for T ~ 72. 5A' while the inset is for T ~ 14.5 
K. Two different scattering rates are used, namely the clean case F = (top curves) and 
F = a/2 R:! 22 meV (bottom curves). For the clean case the curves start at value unity at zero 
frequency. The introduction of impurity scattering reduces the value of the density of states in 
the low frequency range shown in fig. ^ but does not significantly smear the predicted phonon 
structure. Three drops are clearly seen at the frequencies of each of the three delta functions in 
the model spectral density a^F(w) of eq. We note that increasing the impurity scattering 
does not smear out the predicted phonon structure, rather it reduces the effective step drop. 
On the other hand, comparison of the inset with the main frame shows that temperature does 
smooth out the curves considerably, but even for the worst case T = 72.5 K and F = 22 meV, 
the structure in Is{i-ij) remains. For the chosen parameters, which are not extreme, instrument 
resolution does not have a strong smearing effect. We also used a representative extended 
a^F(ijj) with the same A and position of the peaks as in eq. I@J. In this case the results 
(dashed curve) for N{u!) and Is{i^) are indistinguishable on the scale of fig. 

Another way to detect the finite band phonon structures in the EDOS is through tunneling 
spectroscopy. The well known formula for the current, / (assuming the tunneling matrix 
element is independent of energy), through a tunnel junction reads: 
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where V is the applied voltage, e is the charge of the electron and the labels L and R denote 
materials on opposite sides of the junction. The normalization of Ni^^ii[u!) is to its value at 
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Fig. 1 - Comparison of the renormalized density of states N{iL!) (dotted) with the symmetrized 
photoemission signal Is{uj) (soHd) that includes the instrument resolution for a'^F{Lo) of eq. The 
impurity scattering rate is F = and 22 meV for the upper and lower pairs of curves respectively. 
Frequency uj is given in units of half bare bandwidth W/2 = 1.25 eV. The temperature is T = 
72.5 K. For a representative extended a^_F(tj) the results (dashed curve) for N{uj) and Is{uj) are 
indistinguishable on the scale of the plot. Inset: the same as the main frame but for T — 14.5 K. 

Fig. 2 - Differential I-V chacteristic for tunelling between two identical narrow band (i.e. W — 2.5 
eV) metals (lower two curves) and between narrow and wide band metals (upper two curves). The 
temperature is T = 14.5 K (solid curves) and 72.5 K (dashed curves), and the impurity scattering 
rate is F = 22 meV. 



a; = 0. In fig. |2]we show dI{V)/dV vs e *V for a narrow band - wide band case (upper pair 
of curves) and for the case where the same narrow band metal is on both sides (lower pair of 
curves). The impurity scattering rate is F = a/2 « 22 meV and the two temperatures shown 
are T = 14.5 K and 72.5 K. In both cases the phonon structures are clearly identifiable in the 
differential conductance even at higher temperatures. The upper curves are almost identical 
to normalized EDOS N{lu)/N{0) (see inset in fig.^. Features in the lower curves are about 
twice as large in magnitude as in the upper curves. Finite band effects are enhanced when 
the same narrow band metal is used on each side of tunnel junction. 

Optical response. - For characterization of the optical response the quantity of interest 
is the memory function. Its real and imaginary parts are t~j^{uj) and —LuXop{uj) respectively, 
where is the optical scattering rate and Xop is the optical mass renormalization [6]. The 
memory function is connected to the optical conductivity a = ai + ia2 by algebraic relations: 
Tgp^ = {2S/n)U[l/a] and —ujXop — uj ~ (2S'/7r)3[l/f7], where S is the total optical spectral 

weight defined as S' = J^°° dujai{uj). The optical conductivity for our model was obtained 
using linear response theory, neglecting vertex corrections. The result for the real part reads: 

where v"^ is the averaged square of the group velocity (see Ref. [24] for details). We assume 

that the system is isotropic and use for w| the expression u| = derived from 

the quadratic dispersion of free electrons with lower band edge at ^ = —W/2. Here D is the 
number of spatial dimensions, m the free electron mass. Since the complex conductivity <t{uj) 
satisfies the Kramers-Kronig relation, the required imaginary part can be obtained as the 
Hilbert transform of the real part. 
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Fig. 3 - Comparison of the quasiparticle scattering rate — 2E2(w) (dashed) with its optical counterpart 
T~p (w) (soHd) for impurity scattering rates F = 0, 22, and 88 meV, from bottom to top. The 
temperature is T = 14.5 K, the bare bandwidth is = 2.5 eV. 

Fig. 4 - Comparison of — 2Ei(lj) (dashed) with u\op{ui) (sohd) for the same values of F as in fig. El 
but now from top to bottom. The temperature is T = 14.5 K, the bare bandwidth is W = 2.5 eV. 

The memory function is the optical counterpart of the self energy. In fig. Owe compare 
T~^{ijj) with —2Yj2{uj) in the limited range of frequencies oj S [0,500] meV for T = 14.5 K. 
In the infinite band case these quantities are equal at = 0. The bottom pair of curves in 
fig. 121 is for the clean case, the middle pair is for V — a/2 k, 22 meV, while the top pair is 
for r = 2a « 88 meV. The effect of increasing the impurity content is to push the curves 
up corresponding to a larger value of the residual scattering. Beyond this effect there are 
changes in phonon structures, but little smearing. The structures are much more visible in 
the quasiparticle scattering rate — 2E2 (steps) than in the corresponding optical quantity T~p , 
which has only small kinks at frequencies related to the three oscillators of the model a'^F{uj) 
given by Eq. Of course, it is well known that optical data is encoded with the information 
on the phonons even in the infinite band case. However, a second derivative of the data is 
required [15], and this is normally not as accurate. 

In fig. 21 we plot u)\op{ijo) (solid curves) and — 2Si((jj) (dashed curves) for the same values 
of r as those for fig. In this case both quantities become smaller as impurity scattering 
increases in the frequency range shown. Note that all curves cross zero, which is a characteristic 
of finite bands. The crossing frequency depends on T for both Lj\op{ij) and —2Tii{uj), with 
the former dependence stronger. Note that in fig. 01 the phonon structures are again much 
more prominent in the self energy as compared with its optical counterpart. 

Very narrow bands. - So far we have considered the case of a rather wide though finite 
electronic band with full width W — 2.5 eV where W ^ win and the phonon structures in 
EDOS have a generic form. However, band structure calculations for one family of the known 
narrow materials — fulleride compounds — give a much smaller value W — 0.5 eV [25,26]. 
One expects that in this situation many of the discussed features should be significantly 
enhanced [13]. In fig. we present the renormalized EDOS N{lu) as a function of frequency 
w at T = 14.5 K for the model Q!^F(cj) of eq. relevant to KsCqq- The bare EDOS was 
approximated by a flat EDOS with width 0.5 eV, and with fi = (dot-dashed curve). The 
results shown are for impurity scattering rates F — 4.4 and 39.9 meV, as indicated. We used 
two values of the mass enhancement parameter: A = 0.71 [10], which is representative of the 
range quoted in the literature [22,27], and a bigger value A = 2.0, which may be in a range 
more realistic for the fuUerides [28]. 
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Fig. 5 - The renormalized density of states N{lu)/No as a function of frequency ui for the model 
a'^F{u!) of eq. Q relevant to KsCeo- The curves differ by values of A and F as indicated. The 
dot-dashed curve is the bare EDOS with W = 0.5 eV. The temperature is T = 14.5 K. 

Fig. 6 - The optical scattering rate T~p^{uj) as a function of frequency ui for the K^Ceo model of 
eq. ((IJ. Parameters and marking of the curves are the same as in fig. El 



In fig. Owe see that inclusion of the electron-phonon interaction modifies the entire band, 
even for the modest value of A = 0.71. The renormalized EDOS extends well beyond uj = 
W/2 — 250 meV. Note that what determines the energy of the highest steep drop in N{uj) 
is not the bare band edge but rather the energy of the third phonon mode = 190 meV, 
which has the strongest coupling to band electrons in the model of eq. Q). Beyond this 
energy the multiphonon structures appear, with magnitude depending both on A and on F; 
this extends the renormalized band beyond 2.5 x {W/2) . On the other hand, no prominent 
feature corresponding to the bare band edge remains in the renormahzed EDOS. 

The sharp drop in renormalized EDOS in fig.Elat w « is in good quantitative agreement 
with photoemission data on both polycrystalline [27] and monolayer samples [29,30] of metallic 
fuUerides. We note that the model of a^F(w) defined in eq. Q is expected to work reasonably 
well even for two dimensional samples of fuUerides because the stronger coupled W3 mode is 
associated mainly with intramolecular vibrations of Cgo [31]. The flat bare EDOS is also 
appropriate for the two dimensional electonic bands. We conclude that in a very narrow 
band metallic material the electron-phonon interactions have a very profound impact on the 
density of electronic states observable by various spectroscopies, to such an extent that they 
can dominate N{uj). 

The importance of electron-phonon interactions and multiboson processes in understand- 
ing photoemission data in fuUerides has been emphasized by Knupfer et al. [27]. Using a 
different approach they analysed the magnitude of A^(0) and redistibution of electronic states 
derived from the bare band over frequency. In this letter we pointed out the significance of 
finite band effects and provided a useful method to treat other aspects of the problem, such as 
the effect of impurities, within the same formalism. Sensitivity of these effects to the value of 
A could possibly be used to extract the mass enhancement parameter for fuUeride compounds 
from experimental data on the EDOS, supplemented by the other spectroscopic techniques. 
As an example, in fig. we show results for the optical scattering rate in the K^Cqq of eq. Q 
at T = 14.5 K. The curves shown are for A = 0.71 and 2.0, and F = 4.4 meV and 39.9 meV. 
They are marked the same way as in fig. El Note that, similarly to the case of the EDOS, the 
phonon structures dominate the frequency behavior of tZ^{uj) for very narrow band metals. 
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